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ABSTRACT 

Context. The star HD 87643, exhibiting the "B[e] phenomenon", has one of the most extreme infrared excesses for this object class. 
It harbours a large amount of both hot and cold dust, and is surrounded by an extended reflection nebula. 

Aims. One of our major goals was to investigate the presence of a companion in HD87643. In addition, the presence of close dusty 
material was tested through a combination of multi- wavelength high spatial resolution observations. 

Methods. We observed HD 87643 with high spatial resolution techniques, using the near-IR AMBER/VLTI interferometer with 
baselines ranging from 60 m to 130 m and the mid-IR MIDI/VLTI interferometer with baselines ranging from 25 m to 65 m. These 
observations are complemented by NACO/VLT adaptive-optics-corrected images in the K and L-bands, and ESO-2.2m optical Wide- 
Field Imager large-scale images in the B, V and R-bands, 

Results. We report the direct detection of a companion to HD 87643 by means of image synthesis using the AMBER/VLTI instru- 
ment. The presence of the companion is confirmed by the MIDI and NACO data, although with a lower confidence. The companion 
is separated by ~ 34 mas with a roughly north-south orientation. The period must be large (several tens of years) and hence the orbital 
parameters are not determined yet. Binarity with high eccentricity might be the key to interpreting the extreme characteristics of this 
system, namely a dusty circumstellar envelope around the primary, a compact dust nebulosity around the binary system and a complex 
extended nebula suggesting past violent ejections. 

Key words. Techniques: high angular resolution - Techniques: interferometric - Stars: emission-line, Be - Stars: mass-loss - Stars: 
individual (HD 87643) - Stars: circumstellar matter 



Introduction 

Stars with the "B[e] phenomenon" are B-type stars with strong 
Balmer emission lines, numerous permitted Fen lines, and for- 
bidden O i and Fe n lines in their optical spectrum. In addition, 
these stars exhi bit a strong near and mid-IR exces s due to cir- 
cumstellar dust dAllen & Swingslll976r.lContill997l) . For brevity, 
we henceforth refer to stars showing the B[e] phenomenon as 
"B[e] stars". The class of B[e] stars is composed of several sub- 
classes, including objects at different evolutionary stages and of 
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low, medium, and high mass dLamers et al.f l998). Distances to- 
ward them are usually poorly known. As a consequence, the de- 
termination of the evolutionary state of B[e] stars is often uncer- 
tain. 

Observations suggest that the wind of some of the most mas- 
sive Bfej stars (the so-called supergiant B[e] stars or sgB[e] in 
iLamers et alj Il998|) is compo s ed of two distinct components, 
as proposed by IZickgraf et al.l (119851) . The first component is a 
wind of low density in the polar regions; the second component 
is a wind of high density and low velocity located in the equato- 
rial region of the star. The thermal infrared excess is supposed to 
be produced in the outer parts of the equatorial wind where the 
temperature allows the formation of dust grains, but dust might 
also survive much closer to the star in a dense and compact disc. 

The formation of dusty environments around B[e] evolved 
objects is a challenge, with many theoretical issues (see 
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iKraus & La mers 20Ql>. A non-spherical, disc-like environment 
may be the key to understanding dust formation in sgB[e] stars. 
It could be attributed to the rapid rotation of these stars (e.g. 
Lamers & Pauldrachlll99"TtlBiorkman & Cassinelh1ll993tlKrausl 
2006: IZsargo et al.| [2008) but could also be caused by compan- 



ions, since some B[e] stars were found to be binaries. Indeed, the 
number o f companions detected around thes e objects is steadily 
growing ( Miroshnichenko et alj |2004. 2006) , but claims th at all 
B[e] stars are binaries are still controversial dZickgrafl2003l) . 

The number of B[e] stars is rather li mited, and many of them 
are unclassified (cited as unclB[e] in ILamers et al.l Il998l) be- 
cause they exhibit properties usually associated with young and 
evolved objects simultaneously. 

Our group has undertaken a large observing campaign 
to investigate these poorly studied objects (sgB[e] and un- 
clB[e]) mostly using optical interferometry wit h the Very Large 
Telescope Interferometer (hereafter VLTI, see IStee et al.ll2005l 
for an introduction to interferometric observations of Be and 
B[e] star s). Early results were published on the sgB[e] star CPD- 
57 2874 (IDomiciano de Souza et aT1l2007l) . 

HD 87643 (Hen 3-365, MWC 198, IRAS 10028-5825) ap- 
pears to be a s pecial case among sgB [e] and unclB[e] stars. It 
is a B2[e] star dOudmaiier et alj [l998) that exhibits the largest 
infrared excess among this class and appears to be embedded 
in a complex nebula w hose properties are reminiscent of the 
nebulae around LBVs Ivan den Berghll 19721: ISurdei etalJll98U 
I Surdei & Swingslll983l) . HD 87643 is considered a unclBfe] by 
Lame rset al.l (1998): it was classified as sgB [e] (see the esti- 
mation of the bol ometric luminosity in McGregor et al. 1988; 
ICidale et al.l l2001). und er the assumption that i t lies close to the 
Carina arm (2-2.5 kpc [O udmaiiei ^et al.l 1 19981: iMiroshnic henkol 
U998HMcGrego r et al.ll988l). but one can also find it classified as 
a Herbig star dValenti et alj|2000t iBaines et al1l2006l) . The v sin i 
of the central star is unknown since no pure photospheric line is 
detected in the visible spectrum. Clues to the disc-like geometry 
of the environment are provided by the polari sation ellipse in the 
U-Q plane across Ha dOudmaiier et alJll998l). The detected PA 
is about 15°-20° dYudin & Evanslll998t lOudmaiier et al.lfl998l 
On the other hand, Bai nes et al. I d2006h " claimed to have de- 
tected an asymmetric outflow around HD87643 using spectro- 
astrometry. However, they also note that "HD 87643 stands out 
in the complexity of its spectro-astrometry", compared with the 
numerous other Herbig Be stars they observed. 

HD 87643 has displayed a long-term decline of the vi- 
sual brightness for the last 30yrs, from V » 8.5/8.8 in 
1980 (MiroshnichenkoJ[l998), to V ~ 9.3/9.4 in early 2009 
(Poimariski 2009). This decline is superimposed on shorter-term 
variations, such as a 0.5 mag decrease in «1 month, followed by 
a 0.7 m ag increase in »5 months, seen in the ASAS photometric 
survey dPoimanskil2009l) between JD 2453000 and JD 2453 180. 
lYudin & Evansl d 1998b also noted that HD87643 shows Algol- 
like variability, the star showing variations "on a timescale of 
days", and being "bluer during brightness minima". 

The totality of the ISO/SWS 2-45 ytzm mid-IR emission from 
HD 87643 seems concentrated within the smallest ISO aperture 
(14" x 22"), and at lO^m , the dusty environment is unresolved 
on a 1" scale dVoorsI 19991) . The ISO/LWS spectrum also contin- 
ues from the shorter-wavelength SED without any jump, which 
is an additional sign that the mid-IR emission is compact. The 
disc-like geometry is strongly supported by the huge infrared 
excess exhibited by the source, the larg e absorption o f the cen- 
tral star flux, and the polarimetric data. IVoorsI JT999) proposed 
that the disc might be a circumbinary disc. 



In this paper we report new observations that bring a new in- 
sight into this interesting object, proving the presence of a com- 
panion, a resolved circum-primary dust envelope (most likely a 
disc), and circumbinary material. 

The outline of the article is as follows: the observations and 
data recorded are presented in Sect. Q] then we present and dis- 
cuss the main results from our observing campaign in Sect. [2] 
and give a global view of the system, considering this new infor- 
mation, in Sect. [3] 



1. Observations and data processing 

1.1. AMBER/VLTI near-IR interferometry 

HD 87643 was observed at the ESO/Paranal observatory with 
the Astronomical Multi BEam R ecombiner (AMB ER), the near- 
infrared instrument of the VLTI dPetrov et alj|2007l) . The obser- 
vations were carried out on February 18, 2006 in medium spec- 
tral resolution (R = 1500) and during a series of nights in March 
2008 at low spectral resolution (R = 35). AMBER uses three 
8m telescopes (Unit Telescopes, hereafter UT) or three 1.8m 
telescopes (Auxiliary Telescopes, hereafter AT). The calibration 
stars used were HD 109787, HD 86440, HD 101531, HD 63744 
and e Oph. Details of the observations can be found in Table Q] 

Table 1. AMBER and MIDI observing logs. 

projected baseline 



Date 


Stations 


Length 


PA 






[meter] 


[degrees] 


AMBER 1 (UT) 


18/02/2006, 3h 


UT 1-3-4 


95, 57, 129 


22, 88, 46 


18/02/2006, 8h 


UT 1-3-4 


71, 62, 102 


65, 145, 102 


AMBER 2 (AT) 


01/03/2008, 2h 


K0-G1-A0 


83,78, 128 


-168, -94, -132 


05/03/2008, 3h 


G1-D0-H0 


63, 64, 60.0 


-59, 64, 3 


06/03/2008, 5h 


G1-D0-H0 


69, 60, 58 


-37, 90, 19 


06/03/2008, 7h 


G1-D0-H0 


71,52, 33 


-15, 117, 33 


10/03/2008, 3h 


H0-G0-E0 


31, 15.6, 46.8 


-102 


11/03/2008, Oh 


H0-G0-E0 


32, 16, 47 


-147 


12/03/2008, Oh 


H0-G0-E0 


32, 16, 48 


-145 


MIDI (AT) 


26/02/2006, 4h 


D0-G0 


31.3 


76.3 


27/02/2006, 6h 


A0-G0 


57.1 


101.3 


01/03/2006, lh 


A0-G0 


63.5 


38.1 


01/03/2006, 5h 


A0-G0 


61.1 


85.1 


19/04/2006, 4h 


D0-G0 


25.4 


123.1 


23/05/2006, 2h 


A0-G0 


51.7 


120.0 


25/05/2006, lh 


A0-G0 


54.4 


110.6 



The data were processed with the standard AMBER data 
reduc t ion software (amd lib 2.1, see for instance Tatull iet~aT1 
120071: iMillour et a .112004 plus a seri es of advanced scripts to 
calibrate the data dMillour et alj|2008l) . The AMBER DRS per- 
forms a fringe fitting instead of Fourier transforms and computes 
interferometric data products such as V 2 and closu re phases (for 
a review of the interferometri c data analysis, see: Haniff||2(j()ei 
iMonnieri 12006: Millour 2008). We performed the reduction us- 
ing standard select ion criteria (see for instance appendix C in 
Millour et al. 2007) for the individual exposures, rejecting 80% 
of the data before averaging the data products. The additional 
scripts allowed us to compute realistic error bars, including the 
uncertainties on the diameters of the calibration stars, the instru- 
ment atmosphere transfer function instabilities, and the funda- 
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mental noise. The study of the transfer function provided by dif- 
ferent calibrators gives an estimate of absolute errors on V 2 be- 
tween 0.05 and 0. 1 . The visibilities and closure phases are shown 
in Fig. [3] 

1.2. MIDI/VLTI mid-IR interferometry 

The observations of HD 87643 at the ESO/Paranal obser- 
vatory using the MID-I nfrared instrument (hereafter MIDI, 
iLeinert et all I2003L 2004) were carried out from February un- 
til May 2006. MIDI is the mid-infrared (N-band, 7.5-13.5jum) 
two-telescope combiner of the VLTI, operating like a classical 
Michelson interferometer. In our case, only the ATs were used 
since HD 87643 is bright in the mid-IR (156 Jy). W e used a stan- 
dard M IDI observing sequence, as described by Ratzk aet al.l 
(l2007h . The source was observed in the so-called High-Sens 
mode, implying that the photometry from each individual tele- 
scope is performed subsequently to the recording of the fringes. 
The low spectral resolution (R=30) provided by the prism was 
used. The data consist of 7 visibility spectra and one flux spec- 
trum (PSF 1.2" with ATs). The errors on the visibilities, includ- 
ing the internal ones and those from the calibrator diameter un- 
certainty, range from 0.05 to 0.15. The MIDI spectrum was diffi- 
cult to calibrate, and the accuracy of the absolute photometry is 
not better than 30%, although one can see in Fig.BJfhat th e MIDI 
spectrum agrees well with the IRAS dOlnon et alJI 19861) and the 
ISO spectra (affected by probl ems of 'gluin g' between different 
spectral bands, as reported in IVoorsl[l9 99). The log of the ob- 
servations is given in Table Q] We used two different MIDI data 
reduction packages: MIA developed at the Max-Planck-Institut 
fur Astronomie and EWS developed at the Leiden Observatory 
(MIA+EWS0, ver. 1 .5. 1). The resulting visibilities can be seen in 
Fig. a 

1.3. NACO/VLT adaptive optics assisted imaging 

We observed HD 87643 at the ESO/Paranal observatory with 
the NACO adaptive optics camera (N Aos adaptive optics sys- 
tem combined with the COnica camera, iRousset et al.ll2003l) at- 
tached to UT4 of the Very Large Telescope (VLT). NACO was 
operated in the visual wavefront sensor configuration. We ob- 
served the target with Ks (2.2/im) and L' (3.8/mi) broad-band 
filters. The star HD 296986 was used to derive the point-spread 
function (PSF). The S13 camera mode was used for Ks, with a 
13 mas per pixel scale and a 14"xl4" field of view. In L', using 
camera mode L27, the field of view was 28"x28" and the pixel 
scale was 27.1 mas. The auto-jitter mode was used, which, at 
each exposure, moves the telescope in a random pattern within a 
box of side 7" in Ks and 15" in L'. The journal of observations 
can be found in Table|2] 

The data reduction was performed using our own scripts. 
First, bad pixels were removed (i.e. interpolated using the ad- 
jacent pixel values) and a flat-field correction was applied to the 
data. Then, the sky was computed as the median of all expo- 
sures and subtracted exposure by exposure. A visual inspection 
of each exposure allowed us to check the PSF quality. The "bad" 
ones were removed for the next step. This allowed us to signifi- 
cantly improve the final data product image quality compared to 
the standard pipeline-reduced frames. This also led us to disre- 
gard the 20/03/2008 data-sets, which were in any case flagged as 
"failed" in the observing log. A cross-correlation technique was 

1 Available at 
http : / /www . strw . leidenuniv . nl/~nevec/MIDI 



Table 2. Journal of observations with NACO/VLT. All stars 
were observed with a neutral density filter. 



Date 


Mar 


Used 


Used 


Exp 


Seeing 






filt. 


cam. 


time 




zU/U.5/2UU?5, OnzZ 


on Q"7A/1 T 


K-s 


Q 1 1 
ol J 




1 1 T' 

1. Jo 


20/03/2008, 6h36 


HD296986 


K, 


S13 


96s 


1.10" 


20/03/2008, 6h55 


HD 87643 


L' 


L27 


280s* 


1.48" 


20/03/2008, 7h21 


HD296986 


L' 


L27 


324s 


1.03" 


21/03/2008, 3h33 


HD 87643 


K, 


S13 


96s 


0.93" 


21/03/2008, 3h37 


HD 87643 


K, 


S13 


96s 


0.90" 


21/03/2008, 3h43 


HD 87643 


K s 


S13 


84s 


0.91" 


21/03/2008, 3h57 


HD296986 


K, 


S13 


96s 


0.63" 


21/03/2008, 4h21 


HD 87643 


L' 


L27 


315s 


1.14" 


21/03/2008, 4h44 


HD296986 


L' 


L27 


270s 


1.73" 



* saturated 



then used to re-centre the images with about 1 pixel accuracy. 
Finally, all the selected frames were co-added, resulting in the 
total exposure time shown in Table [2] 

As a result, we got a pair of science star/PSF star images for 
each observation (which were repeated due to changing weather 
conditions and saturation of the detector). The PSF FWHM is 
75+4 mas in the K v -band and 111+6 mas in the L'-band. In addi- 
tion to a deconvolution attempt presented in Sect. 12.21 we com- 
puted radial profiles to increase the dynamic range from about 
» 10 3 per pixel to * 10 4 -10 5 . The result can be found in Fig. [5] 

1.4. 2.2m/WFI large-field imaging 

We retrieved unpublished archival ESO/Wide-Field Imager 
(WFI) observations of the nebula around HD 87643 from the 
ESO database^ carried out on March 15 and 16, 2001. The WFI 
is a mosaic camera attached to the ESO 2.2m telescope at the 
La Silla observatory. It consists of eight 2kx4k CCDs, forming 
an 8kx8k array with a pixel scale of 0.238" per pixel. Hence, 
a single WFI pointing covers a sky area of about 30' x30' . The 
observations were performed in the B, V, R c and I c broadband 
filters and in the narrow-band filter in the Ha line (A = 658nm, 
FWHM = 7.4nm). However, we used only the observations 
with B, V, and R fi lters given that it is a reflection nebula only 
(ISurdei et al.ll 1 98 lb - In order to cover the gaps between the WFI 
CCDs and to correct for moving objects and cosmic ray hits, for 
each pointing, a sequence of five offset exposures was performed 
in each filter. 

The data reduction was performed using the package 
alambic developed by B. Vandame based on tools available 
from the multi-res olution visual model package (MVM) by 
iRue & Biiaouil d!997l) . The image is shown in Fig. [6] 

2. New facts about HD 87643 

2. 1 . interferometry data: A binary star plus a compact dusty 
disc 

From the sparse 2006 AMBER medium spectral resolution data, 
the structure of the source could not be inferred unambiguously. 
The squared visibilities from 2006 did not show any significant 
variation with spatial frequency. The closure phases were mea- 
sured as zero within the error estimates (see Fig. [3j. 

2 http://archive.eso.org 
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The Bry emission line appears clearly in the AMBER 
calibrated spectrum after using the technique described in 
lHanson et al.l (l996) of fitting and removing the Bry line for the 
calibrator. The resulting Bry line accounts for 15% of the con- 
tinuum flux at its maximum and is spread over 5 spectral pixels 
250k m • s _1 ), similar to what is reported in McGregor et al. 



(250 
(1988). 



1 « / 
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% 



✓ 



- sp. freq. (cycles/as) 



-200 

- sp. freq. (cycles/as) 



Fig. 1. UV coverage obtained on HD87643 in the H-band (left), 
in the K-band (right). The radial lines are caused by the different 
wavelengths within the bands. 



The extensive AMBER low spectral resolution data, secured 
in 2008 (see Fig. |3J, show large variations of V 2 and closure 
phases and clearly indicate a modulation from a binary source. 
Given the good uv coverage (see Fig. [T]i and the quality of the 
data, we undertook an image reconstruction process, aiming to 
determine the best models to use for the interpretation. 

2.1.1. Image reconstruction 

The images in Fig. [2] were reconstructed from the AMBER 
squared visibili ties and closure phases using the MIRA software 
(Thiebaut 2008). MIRA compares the visibilities and the clo- 
sure phases from a modelled image with the observed data using 
a cost-estimate optimisation, including a priori information such 
as image positivity (all image pixels are positive) and compact- 
ness of the sou rce (using a so-called L2-L1 regularisation, see 
lThiebautll2008L for details). 

Separate image reconstructions were made using the H-band 
and K-band data, assuming that the object was achromatic over 
the whole band considered. MIRA needs a starting image, but we 
found that the image reconstruction result did not significantly 
depen d on it. Extensive tests and comparisons with oth er soft- 
ware dHofmann & Weigelt|[T99ll: iBaron & Youngll2008L see ap- 
pendix[A} made us confident that we could distinguish between 
artifacts and real structures in the images. The theoretical field 
of view, computed using the formula 2.44/l/B m j n , is ss 70 mas, 
B m i n being the shortest projected baseline (15 m). We convolved 
the K-band image with a Gaussian beam of FWHM A/B msix — 
3.5 mas and the H-band one with a beam of A/B max = 2.7 mas, 
B max being the longest projected baseline (128 m). The image 
noise is of the order of 1%. 

The reconstructed image reveals the presence of a com- 
panion star with a plane-of-sky separation of x;34.5mas (see 
Table [3]). the K-band image also reveals an extended structure 
detected around the main star. We cannot assess any elonga- 
tion of this resolved structure due to the poor UV coverage in 
the NW direction. We emphasise that there is still a 180° un- 



certainty due to the overall uncertainty on the AMBER closure 
phase sign (and not due to the image reconstruction process: see 
appendix lA.1.3b . We will refer in the following to "the northern" 
and "the southern" components. 

We also extracted relative fluxes for the different components 
within the dashed circles in Fig. [2] (see Table [3), but one has to 
bear in mind that extracting such fluxes is affected not only by 
large errors, but also potentially by large systematics. Therefore, 
we also performed the flux measurements using model-fitting. 

2.1.2. Model fitting 



In Section l2~m . we show that the 2008 AMBER data can be di- 
rectly interpreted assuming the presence of a binary source with 
a plane-of-sky separation of «34 mas and a RA. of w 0° (see 
Table [3} and a resolved southern component. We fitted the data 
using such a model to obtain more accurate component fluxes. 
The images show a series of large-scale artifacts (i.e. flux whose 
spatial location is poorly constrained), accounting for a large 
fraction of the total integrated flux. These artifacts are related 
to the lack of data with high visibilities at small spatial frequen- 
cies. We therefore also included an extended component to our 
model, fully resolved on all AMBER baselines. 

Our model-fitting used the scientific software yoriclfl com- 
bined with the AMBER data reduction software amdlib. This 
tool was complemente d with a serie s of optimisation scripts de- 
veloped by the JMMC dBechet et al.ll2005l) and others developed 
by us. We judge that the resulting fits were satisfactory, even 
though the formal^ 2 is about 40 (see Fig. [3}, probably due to an 
underestimation of the errors during the data reduction. 




10 20 

Sp. Freq (B pra A, cycles/arcsec) 

Fig. 4. MIDI visibilities (in Gray with error bars), projected on the 
direction of the detected binary star and compared with our model (in 
thick black line). 

The MIDI data show complex, spectrally dependent visi- 
bility variations from one observation to another (see Fig. 0). 
Spherical and even 2-D axi-symmetric dusty models were not 
able to account for this data set. We used the separation and P. A. 
inferred from the AMBER 2008 data as initial values for the fit- 
ting of the AMBER 2006 and MIDI 2006 measurements. 

The result of such a fit shows that our binary model is 
compatible with both the AMBER and MIDI 2006 data sets. 
Moreover, we find that the MIDI separation is close to that 
derived from the AMBER 2008 data set. The P.A. of the 
secondary component differs slightly between near- and mid- 
infrared, which might be due to an offset of the dust structure 



3 open-source scientific software freely available at the following 
URL: http : //yorick . sourcef orge . net 
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Fig. 2. Aperture-synthesis images of HD87643 reconstructed from the H-band AMBER data (panel[2^) and the K-band data (panel[2p), assuming 
the source is achromatic in each band. Contours with 50, 10, and 1% of the maximum flux are shown, and the beam size is shown in the lower-right 
box. The arrow shows the companion star, whereas the dashed circles show the trusted structures in the images. 



2 0.4 - 



> 

T3 



C/2 0.2 




50 100 150 

Sp. Freq (BpraA, cycles/arcsec) 



50 100 150 

Sp. Freq largest base (B pr0 j/X, cycles/arcsec) 





Sp. Freq (Bp ro j/X, cycles/arcsec) 



50 100 150 200 

Sp. Freq largest base (B pr0 ]/X, cycles/arcsec) 



Fig. 3. Squared visibilities (left) and closure phases (right) plotted as a function of spatial frequency, projected onto the inferred binary position 
angle. The 2006 AMBER data set is seen on top, 2008 on the bottom. The model involving a binary system composed of one resolved and one 
unresolved component plus a fully resolved background is shown with a solid black line. 



compared to the main star. The total flux seen by MIDI is dom- 
inated by a fully resolved background in addition to the unre- 
solved binary components (TableO. The AMBER 2006 data set 
is too limited to provide a strong constraint on the the flux ratios; 
therefore, we do not use them in the following analysis. 



We stress that the MIDI data is affected by an independent 
180° orientation ambiguity to the AMBER ambiguity (due to the 
uncalibrated closure phase sign), as it only consists of squared 
visibility data. Hence, we cannot definitely relate the compo- 
nent fluxes derived from the MIDI data to those derived from 
the AMBER data. 
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Table 3. Estimated diameter of the envelope, position of the secondary source, and respective fluxes, relative to the total flux. 
The envelope flux includes both the southern component and its compact circumstellar envelope fluxes. The errors given for the 
model-fitting reflects only formal estimates from the x 2 value. 





AMBER (model fitting) 


AMBER (Image reconstruction) 


MIDI (2006) 


Wavelength (yum) 


1.55 I 2.45 


H-band | K-band 


8.1 


10.5 | 12.9 




Point + Gaussian + Background 


Image 


Point + Gaussian + Background 


Envelope size (mas) 


3.1+0.1 


4.5+0.1 






< 6.5 


< 10.6 


< 3.5 


Envelope flux (% total flux) 


61+2 


72 + 2 


ss 46 


-59 


17 + 2 


9+1 


6.3+0.5 


Background flux (% total flux) 


12 ±2 


9 + 2 


* 44 


» 25 


69+ 13 


77 + 9 


75 + 9 


Secondary flux (% total flux) 


27 ± 1 


19+1 


» 10 


» 16 


14 + 3 


15 + 2 


19 + 2 


Secondary separation (mas) 


33.48 ± 0.03 


34.5 ± 0.5 


34.6 ± 0.5 


36.8 + 0.1 


Secondary PA (°) 


-4.50 ± 0.03 


-2.5+0.8 


-4.1+0.8 


-26.1 ± 1.7 



2.2. NACO K-band imaging: a very compact dusty 
environment 

We performed a deconvolution process on our NACO images. 
Using both speckle techniques and the Lucy-Richardson decon- 
volution algorithm, an elongated north-south compact structure 
could be detected in the K-band but not in the L-band. This is 
compatible with the AMBER data. 

In addition, we tested the presence of dusty emission in the 
surroundings of the star by azimuthally integrating the flux to 
increase the dynamic range (see Fig. |5j. In the K-band, no ex- 
tended emission is detectable up to 1.5 arc-seconds (at larger 
distances, ghost reflections or electronic ghosts impair the dy- 
namic range of the images) and a faint emission can be seen in 
the L-band from 0.5 to 3.0 arc-seconds. This might be linked 
to the fully resolved component seen in our AMBER data, and 
therefore probably implies that it is much more extended than 
the field of view of AMBER. However, the dynamic range of 
our images (« 10 3 per pixel, and a 10 4 -10 5 for the radial pro- 
files) does not allow us to access the spatial distribution of such 
a dusty nebula. 

2.3. WFI imaging: The large scale nebula revisited 

The extended nebu la was disc o vered and studied by 
van d en Berghl d 19721) : ISurdei et all dl98ll) : ISurdei & Swings! 
dl983). In particular, using the 3.6m ESO telescope at La 
Silla, Surdei & Swingsl d 19831) showed that the structure was a 
reflection nebula with a 60-100" extension. 

The nebula revealed by the WFI instrument (see Fig. [6} with 
the 2.2m telescope at La Silla shows increased dynamic range 
but shares the same spatial resolution and global morphology 
as the previous works. The nebular structures are primarily seen 
in the north-west quadrant from the central star in the form of 
an extended and structured filamentary feature at distances of 
15-20" (north) to 35-50"(west). In some places, labelled (Bl) to 
(B4), the nebula appears blown-up by the central star wind, with 
clumpy and patchy features. 

This structure has no extended counterpart in the south-east 
region, where arc-like structures (labeled (Al) to (A4) in the fig- 
ure) can be seen closer to the star (at 10", 12", 14", and 18" 
respectively). These arcs are spaced by two to four arc-sees and 
might come from past ejection events. 

In the south-west region, at 12", 19", and 30" respectively 
from the star, two bright and one faint arc can be seen (labelled 
(A5) to (A7) in the figure). These arcs do not appear to be con- 
nected to the previously-mentioned ones. 

A star count made using sextractor dBertin & Arnoutsl [l996l) 
gives values twice as low in the east half than in the west half 
of the image. No significant relative change in the counts as a 



function of filter is detectable over the field. This indicates that 
there is a total absorption screen somewhere in the field, totally 
masking the background stars in the east half of the WFI image. 
We note that it corresponds to a dark cloud in the eastern part, 
clearly seen in our WFI image in front of red (hydrogen ?) neb- 
ular emission. However, we cannot directly assess if this screen 
is in front or behind HD 87643 and its nebula. 



3. Discussion 

3. 1 . Adopting a distance for the system 

HD 87643 is thought to be an evolve d Bfel star, and the neb - 
ula suggests a link to the LBV class (ISurdei & Swingsl [l983). 
HD 87643 stands out as unusual compared to typical sgB [e] stars 
such as CPD-57 2874, since ISO SWS and LWS spectra show a 
considerable amount of cold dust (T<150K). 

iMcGregor et all dl988l) implicitly assumed that HD 87643 
is a supergiant in order to determine its distance, based on its 
location in the direction of the Carina Arm. Such a location 
would place it at a distance of 2-3 kpc. However, the distance 
of HD 87643 has never been accurately measured. When other 
estimates of the distance are found in the literature, the values 
are usually much smaller. These were derived by: 



1. Ivan den Berghl d 19721) and Ide Freitas Pacheco et all d!982l) : 
530 pc, based on Kurucz models for a B2V star (T e jf = 
18000 K, E(B-V) = 0.63 mag, R« = 6Rq); 

2. ISurdei et ail d!981l) : 1 kpc, based on M v = -4. 1 derived from 
photometry ; 



Shore et al. ( 1990): 1.2 kpc, using the IUE UV spectrum; 



3. 

4. lLopes et al.l(ll992l) : 2-3 kpc using the equivalent width (here- 
after, EW) of Nai lines. We note, however, that the high 
reported value of 2.9 kpc is affected by a factor of 2 error 
in thei r application of the statistical reiation cited by lAllenl 
(11973b . Considering the correct relation (r = 2.0D, r be- 
ing the distance in kpc and D being the mean EW in A of 
the two Na D lines), the distance would be 1 .46 kpc (Lopes, 
private communication). The EW were measured using low- 
resolution spectra, implying a possible contamination of cir- 
cumstellar origin ; 

5. lOudmaiier et al.l (Il998l) : 1-6 kpc, also based on the Nai D- 
line equivalent widths derived from a higher-resolution spec- 
trum, but using data on the interstell ar extinction for nearby 
stars instead of the relationship from Allen] d 1973b . The large 
uncertainty of this distance is due to the large scatter of the 
extinc t ion va lues in the Carina arm line of sight. 

6. IZored (fT998): 1 .45 kpc, based on a SED fitting using various 
assumptions. 
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Fig. 5. NACO radial profiles. Left is the K, filter and right is the L filter. HD 87643 is shown as full lines (3 observations in Ks, and one in L', 
see Table[2}. The calibrator star is shown as a dotted line. 




a (as) alpha(as) 

Fig. 6. Reflection nebula around HD 87643, as a composite of R, V and B filters (left), together with a sketch presenting the main structures (right). 
The saturated regions are masked by black zones in the image. This image is a small part of the 30'x30' field-of-view of the WFI observations. 
The sketch shows HD 87643 as a red circle, whereas other stars are marked as red crosses. The nebular contours are drawn as black lines and the 
prominent features are labelled (Al) to (A7) for arc-like structures and (Bl) to (B4) for apparently blown-up nebular structures. Structures which 
are faint or uncertain are marked as dotted lines. 



Therefore, we arbitrarily use the distance of 1 .5 kpc, as sug- 
gested by some of these works, in the subsequent parts of the 
present article. 

3.2. Orbital characteristics of the binary system 

One of the main results from the AMBER interferometric obser- 
vations is that HD 87643 is a double system. In the K-band, the 
NACO image core is also significantly resolved with an elonga- 
tion in the same direction as the AMBER binary. The separation 
and position angle of the components are well constrained by the 
data. It is more difficult to derive an accurate flux ratio, given the 
complexity of the object. Here, we shall summarise our results, 
from the most constrained ones to the most speculative. 

- The plane-of-sky separation of the components is 34 ± 
0.5 mas. At 1.5 kpc, it corresponds to a projected separa- 
tion of 51 ± 0.8 AU. As a comparis on, the separation of the 
interacting components of f3 Lyrae dZhao et al.l 12008) is of 



the order of 0.2 AU (orbital period 12 d), the separation of 
the y 2 Velorum (Millour et al. 2007) colliding-wind compo- 
nents is of the order of 1 .3 AU (period 78 d), and t he approxi- 
mate s eparation of the high - mass binary 6 l Ori C dStahl et al.l 
120081: iPatience et al.l 120081: iKraus et all 120071) is 15-20AU 
(orbital period of between 20 and 30yrs). 

- Therefore, the orbital period of HD87643 is large, approxi- 
mately 20-50 years. In particular, no significant orbital mo- 
tion is seen between our 2006 and 2008 data sets. 

- The orbital plane might be seen at a high inclination (i.e. 
close to edge-on). Indeed, thi s is suggested by the high 
level of observed polarisation (Oudm aiier et al.l ll998). The 
large-scale measurements also seem to imply a bipolar 
morphology (see Sect. |3.3I >. In addition, HD87643 shows 
photometric variability. The short-term variations (e.g. am- 
plitude ~0.5mag, seen in the ASAS light curves and in 
iMiroshnichenkoll 1 9981) are similar to Algol-like variability 
dYudin & Evanslll998l) and are therefore probably a conse- 
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quence of the time-variable absorption by the material pass- 
ing through the line-of-sight. This is also in favour of a 
highly inclined (edge-on) disc-like structure. 
- The orbit might be highly eccentric. This hypothesis is sup- 
ported by the periodic structures seen in the large-scale neb- 
ula (see Sect. 13.3b that would suggest periodic eruptions. 
These eruptions would be a sign of close encounters of the 
binary components, while we observed a well-separated bi- 
nary. 

3.3. A link to the larger-scale nebula 

Dust is found at large distances east and west of the nebula, with 
fluxes reaching 15 Jy at 16.7 /mi (ISO/CAM1, offset of 25" from 
the star, aperture 20"xl4"), whereas no flux is detected in the 
north and south position^. This points to a bipolar nebulosity in 
the east-west direction. 

As mentioned by Surd ej et all d 198 lh . the number of field 
stars decreases in the south-east direction. Also from this work 
(Sect. 12.31 ), the number of visible stars in the eastern part of the 
WFI image of the nebula is about 2 times less than in the western 
part. This may mean that the other side of a symmetrical bipolar 
nebula remains hidden in our WFI image. 

Finally, we note that inferring the geometry of the nebula 
is potentially of great importance as this reflection nebula al- 
lows one to study the central star from different viewing angles, 
i.e probing different latitudes of the ce ntral star. The aniso tropy 
of the star flux was alr eady noted in ISurdej et al.1 (1 198 lh and 
ISurdei&Swingsl(fT983h . 

lOudmaijer et alJ (Il998h and iBaines et all (120061) measured 
expansion velocities of ss 1000 km • s _1 . Therefore one can es- 
timate an ejection time for the nebula: at 1.5 kpc distance and 
taking a 50" extent, one gets »355 yrs. 

For the arc-like structures and the same value for the expan- 
sion speed and distance, we find the following values: «71 yrs 
for Al, *85 yrs for A2, « 100 yrs for A3, and «128 yrs for A4. 
This first series gives ejection time intervals of a; 14, «14 and 
a;28 yrs, respectively, between two consecutive arcs. This might 
be the trace of a periodic ejection with a period of a; 14 yrs on 
the assumption that one arc (between A3 and A4) is not seen 
in our image. However, this periodicity is probably affected by 
a projection effect since the vsinz of the arcs is not known and 
they appear to be almost linear (instead of circular); hence, the 
ail 4 yrs periodicity is a lower limit. Concerning the second series 
of arcs, we find «85yrs for A5, ~135 yrs for A6, and «213yrs 
for A7. In this case, the apparent ejection time intervals between 
two consecutive arcs are w50yrs and »78 yrs. Given that the arcs 
appear almost circular in the image, we can assume that the pro- 
jection angle is close to 90°. Since we have only three arcs in 
this case (with one barely seen), we can only put an upper limit 
of ^50 yrs on the periodicity of these ejections. 

These broken structures suggest short, localised ejection that 
might coincide with short periastron passages of the compan- 
ion, triggering violent mass-transfer between the components. 
Therefore, we tentatively infer limits between 14 and 50 yrs for 
the periodicity of the binary system. Monitoring the system at 
high angular resolution over a timescale of a few dozen years 
would most likely bring an unprecedented insight into this sys- 
tem. 



4 These measurements can be found in the ISO database: 
http : //iso . esac . esa . int 



3.4. The nature of HD 87643 

Our interferometric measurements show a complex object com- 
posed of a partially resolved primary component, a compact 
secondary component, and a fully resolved component (i.e. ex- 
tended, or nebular, emission). As shown in Table[3] their relative 
flux strongly varies between 1.6 and 13jum. 

The L-band NACO image of HD 87643 cannot be distin- 
guished from the calibrator star, and there is very little emis- 
sion at l"-5" distance. In the K-band, except from the binary 
signature, no flux can be detec ted at a larger distance . Given the 
amount of dust in the system dMcGregor et alj[l988l) . the com- 
pactness of the near-IR emissions indicate that most of it resides 
in a region smaller than aTOOmas. We may assume that the K- 
band flux comes entirely from the very central source, as seen by 
NACO (i.e., it contains the central binary star plus the extended 
emission, as seen by AMBER). We assume the same for the H- 
band. Given these hypotheses, we can infer the absolute flux of 
each component using (for example) the 2MASS magnitudes in 
the H and K-bands. 

In the mid-infrared, the MIDI spectrum is close to the ISO 
one. Therefore, we can also assume that all the N-band flux orig- 
inates from the MIDI field of view (i.e ~1") and infer the abso- 
lute fluxes of the components in the N-band. These fluxes (N- 
band) and magnitudes (H and K-bands) are presented in Table [4] 
and plotted in Fig.|7]as a function of the wavelength. 



Table 4. Estimated fluxes for each component from the model 
fitting of Sect. 12.1.21 We used 2MASS magnitudes as the total 
magnitude in H and K-band and the MIDI fluxes in N-band. 



Band or A Primary Background Secondary 

(total value) (W/m 2 /jjm) (W/ in 2 / fim) (W/m 2 /jim) 



H (4.8 mag) 


8.5 x 10- 


12 


1.6 x 10- 


12 


3.7 x 10- 


-12 


K (3.5 mag) 


1.2 x 10- 


-ii 


1.5 x 10- 


-12 


3.1 x 10- 


-12 


8.1/«n (92 Jy) 


7.2 x 10~ 


-13 


2.9 x 10- 


■12 


5.9 x 10- 


■13 


10.5yum(146Jy) 


3.5 x 10- 


-13 


3.0 x 10- 


■12 


6.0 x 10- 


-13 


12.9yum(126Jy) 


1.4 x 10- 


-13 


1.7 x 10- 


12 


4.3 x 10- 


-12 



Even if the effect of reddening on the H and K-band flux is 
not negligible, a significant reddening fraction comes from the 
circumstellar envelope. Therefore we performed this study on 
the original measurements, without de-reddening. 

3.4.1 . A dust-enshrouded star in the South 

The flux of the envelope around the southern component, from 
the AMBER and MIDI data, can be qualitatively described by a 
1300K black body radiation (dotted line in Fig.[7]l. Its Gaussian 
FWHM, inferred from the model-fitting, is ~4mas (i.e. 6 AU at 
1 .5 kpc) in the near-infrared (H and K-bands), and is well con- 
strained by the AMBER 2008 data. These suggest an extended 
dusty envelope around this component and clearly indicate that 
the dust must be very close to the sublimation limit. Thus, the 
H and K-band emission should mainly originate from the in- 
ner radius of a dus ty disc , encir cling a viscous gas disc or a 
2-component wind (Porter 2003). Estimating the envelope ex- 
tension with a ring instead of a Gaussian would lead to a slightly 
smaller size. Thus, we can estimate that the inner radius of the 
dusty disc is of the order of 2.5-3 AU (at 1.5 kpc). 
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Fig. 7. A view of the (non-de-reddened) HD87643 SED with 
the extracted fluxes from our interferometric measurements. The 
southern component flux is shown in red (top curve at 2 fim), 
the northern component flux is shown in blue (middle curve at 
2 fim), and the resolved background flux is shown in pink (bot- 
tom curve at 2 fim). The dotted lines correspond respectively 
to a Kurucz spectrum of a B2V star and to a black-body flux at 
1300K, for comparison. 



3.4.2. A puzzling dusty object in the North 

By contrast, the northern component is unresolved (size, or di- 
ameter, of the source < 2 mas, i.e. < 3 AU at 1 .5 kpc). Moreover, 
its flux can be accounted for neither by a simple black-body at a 
constant temperature nor by free-free emission. The slope of the 
flux variation between the H and N-bands suggests a range of 
temperatures for the dust (from at least *300K to »1300K). No 
cold dust can survive so close (Tx;300K at a radius < 1.5 AU) 
to a putative luminous hot star without effective screening of the 
stellar radiation. 



3.4.3. A cold dust circumbinary envelope 

The resolved component shows a large increase of flux between 
H, K, and N and also carries most of the silicate emission (see 
Fig. [7]) . The binary system is therefore embedded within a large 
oxygen-rich, dusty envelope whose shape is not constrained by 
our data. This envelope must significantly contribute to the red- 
dening of both components. 

In conclusion, we propose the following picture of the sys- 
tem: two stars with dusty envelopes probably surrounded by a 
common dusty envelope. One source might be a giant or a super- 
giant hot star surrounded by a disc, in which we would mainly 
see the inner rim corresponding to the dust sublimation radius. 
The other source is surrounded by a compact dusty envelope, 
and it is either not a luminous hot star or it is heavily screened 
by very close circumstellar matter. In the N-band, the circumbi- 
nary envelope contribution is not greater than * 1", since the 
MIDI (aperture 1.2") and ISO (aperture 22x14") fluxes match 
very well, and is not smaller than « 200 mas, as our MIDI data 
indicates. 



4. Concluding remarks 

Our work presents new observations of HD87643, including 
very high angular resolution images. It completely changes the 
global picture of this still puzzling object: 



- The binary nature of the system has been proved, with a pro- 
jected physical separation of approximately 51 AU in 2008 
at the adopted distance of 1 .5 kpc. The orbital period is most 
likely several tens of years, and the large-scale nebula indi- 
cates a possible high eccentricity. 

- The temperature of the southern component is compatible 
with an inner rim of a dusty disc. In the near-IR, we do not 
see the central star, but it might be a giant or a supergiant 
star. 

- The northern component and its dusty environment are un- 
resolved. The underlying star is unlikely to be a massive hot 
star, or is heavily screened by close-by circumstellar mate- 
rial, as cold dust (T*300K) exists closer than 1.5 AU from 
the star (if it lies at 1.5 kpc). 

- The system is embedded in a dense, circumbinary, and dusty 
envelope, larger than 200 mas and smaller than 1". 

High angular resolution observations were definitely vital in 
partly revealing the nature of this highly intriguing object. 

Therefore, we call for future observations of this system, us- 
ing both high spectral resolution spectroscopy and high angular 
resolution techniques, to place this interesting stellar system on 
evolutionary tracks and better understand its nature. 
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Appendix A: MIRA imaging of HD87643: tests and 
reliability estimates 

This work made use of the MIRA software to reconstruct images 
from the sparsely sampled AMBER interferometric data. In this 
appendix, we present test results demonstrating which structures 
in the reconstructed images are reliable. 

A.1. Altering the UV coverage 

The idea here is to find out if the structures remain qualitatively 
the same in the final image when altering the resolution (remov- 
ing the high spatial frequencies) or the low frequency compo- 
nents of the UV plane. 

A.1 .1 . Removing low spatial frequencies 

Removing low spatial frequencies reduces the field of view, and 
the MIRA image reconstruction software does not manage to 
overcome this difficulty (see Fig. IA.ll . Therefore, we conclude 
that the low spatial frequencies (i.e. the short baselines) are just 
as important as longer baselines to perform image reconstruc- 
tion. 




E < sp. freq. (cycles/as) 



Fig.A.l. Removing low spatial frequencies to check whether 
the structures are reconstructed. 



A.1.2. Removing high spatial frequencies 

As a second test of the image reconstruction reliability, we cut 
all the high spatial frequencies in order to get a more symmetric 
UV coverage. The result is seen in Fig. IA.2I together with the 
corresponding UV coverage. The striking point compared to the 
image presented in Fig.|2]is that all the structures seen in this new 
image are qualitatively the same as before. This ensures that the 
binarity and resolution of the southern components are not due 
to image reconstruction artifacts. 




E < sp. freq. (cycles/as) 



Fig.A.2. Removing high spatial frequencies to have a more 
symmetric UV coverage. 



A.1.3. Imaging without phases 

Since the AMBER closure phases in our data-set were noisy 
(<r ss 0.5 - 1 radian), we wondered whether these phases were 
contributing significant information to the image reconstruction 
compared to squared visibilities alone. 

We tried to reconstruct an image of HD87643 using the same 
AMBER data, except for the closure phases, and both the MIRA 
and BSMEM image reconstruction software. Indeed, they are 
able to cope with phase-less data, in this case making a phase- 
retrieval image reconstruction. In Fig IA.3I we show the MIRA 
reconstruction. We are able to reconstruct qualitatively the same 
image of HD87643 as in Fig. [2] without using the closure phase 
information. 

We note, however, that we had to rotate the result to match 
the orientation of Fig. [2] This makes sense since only the phase 
information would be able to orient the system. 
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Fig. A.3. HD87643 image reconstruction without the phase in- 
formation. We had to rotate the image by 180° to make it match 
the previous reconstructions. 

A.2. Testing with other image reconstruction algorithms 

MIRA is only one example of an image reconstruction software 
package for optical interferometry. Other software packages ex- 
ist, and we present here the tests we made using several of them. 

A.2.1. BSMEM 

BSMEM is based on the Maximum Entropy Method (MEM) ap- 
plied to bispectrum measurements (deduced from the squared 
visibilities and closure phase measurements). It works using 
an iterative algorithm, comparing the distance from the recon- 
structed image to the data an d the "entropy" com puted from the 
properties of the image itself (Baron & Young 2008). We applied 
BSMEM to the AMBER data. The result can be seen in Fig. |A.4| 
We find the same structures as for the MIRA reconstruction pro- 
cess. 
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Fig. A.4. Image reconstruction of HD87643 using the BSMEM 
software. 



at the particular coordinate, which leads to a new model image 
that minimises the deviations (x 1 ) between the model bispectrum 
and the measured object bispectrum elements. An approxima- 
tion of the ^-function was derived which allows fast calculation 
of a large number of iteration steps. Adding both positive and 
negative building blocks, taking into account the positivity con- 
straint, and adding more than one building block per iteration 
step improves the resulting reconstruction and the convergence 
of the algorithm. The final image is obtained by convolving the 
Building-Block reconstruction with a beam matching the max- 
imum angular resolution of the interferometer. We applied this 
method to the AMBER data, and the resulting image (Fig. IA.5t 
shows many similarities with the MIRA one, including a re- 
solved southern component in the K-band. 
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Fig. A.5. Image reconstruction of HD87643 using the Building- 
Block software. 



A.3. Conclusion 

The tests presented above show that the binary and resolved 
southern component are structures that can be trusted in the im- 
ages. All other structures (possible elongation of the southern 
component, inclined large "stripes" in the images) are artifacts 
from the image reconstruction process, which may indicate that 
the image has some additional flux, not constrained by the ob- 
servations (fully resolved background). 



A.2.2. Building-Block Method 

The Building-Block method ( Hof mann & Weigeld 19931) was de- 
veloped to reconstruct diffraction-limited images from the bis- 
pectrum of the object obtained with bispectrum speckle inter- 
ferometry and long-baseline interferometry. Since the intensity 
distribution of an object can be described as a sum of many 
small components, the Building-Block algorithm iteratively re- 
constructs images by adding building blocks (e.g. 5-functions). 
The initial model image may simply consist of a single 6 peak. 
Within each iteration step, the next building block is positioned 



